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ABSTRACT 

An investigation of translaminate fracture and a progressive damage methodology was 
conducted to evaluate and develop residual strength prediction capability for laminated 
composites with through penetration notches. This is relevant to the damage tolerance of 
an aircraft fuselage that might suffer an in-flight accident such as an uncontained engine 
failure. An experimental characterization of several composite materials systems revealed 
an R-curve type of behavior. Fractographic examinations led to the postulate that this crack 
growth resistance could be due to fiber bridging, defined here as fractured fibers of one ply 
bridged by intact fibers of an adjacent ply. 

The progressive damage methodology is currently capable of predicting the initiation 
and growth of matrix cracks and fiber fracture. Using two different fiber failure criteria, 
residual strength was predicted for different size panel widths and notch lengths. A ply 
discount fiber failure criterion yielded extremely conservative results while an elastic- 
perfectly plastic fiber failure criterion showed that the fiber bridging concept is valid for 
predicting residual strength for tensile dominated failure loads. Furthermore, the R-curves 
predicted by the model using the elastic-perfectly plastic fiber failure criterion compared 
very well with the experimental R-curves. 
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CHAPTER I 
INTRODUCTION 

Industrial Applications and Problems 

Damage tolerance has recently been the underlying issue driving much of the research 
in aircraft structure design. Consider, for example, any modem civilian aircraft with a 
large number of accumulated flight hours. For such an aircraft the flaws in a particular 
component may become critical before detection. Such a component will have to be 
designed to a safe life because it is not damage tolerant The structural component will 
have to be overhauled or replaced before the aircraft completes a certain number of flight 
hours. A damage tolerant design insures detectable damage before catastrophic failure. 
Damage tolerance is the attribute of a structural component such as the fuselage or wing that 
allows the aircraft to survive an in-flight accident such as an uncontained fan blade or 
engine failure, or an impact with a foreign object such as a bird. 

A familiar incident that occurred in 1988 that escalated damage tolerance research was 
the inflight structural failure of a particular Airlines flight [1,2], Figure 1. The upper 
fuselage ripped open and a large section of the skin peeled away. The cause for this failure 
has been identified as a flaw in the rivet design. The stress concentrations at the knife 
edges of the rivet holes caused the formation of small fatigue cracks. Since these cracks 
occurred between the skins of the lap splice joints, they were virtually undetectable. To 
insure that this didn't occur with any of the other 737's, doublers were inserted using a hot 
bond adhesive between the riveted joints. Activities to insure the safe operation of the 
aging aircraft fleet includes increased maintenance and inspections, repair and 
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modifications, new advanced inspection technology to locate visually undetectable fatigue 
cracks, and research to design damage tolerant materials and structures. 

Damage tolerance studies have discovered similarities between metals and fiber- 
reinforced composites as well as some obvious differences. Unlike homogeneous metals, 
composite materials are made of multiple constituents. Therefore, the modes and types of 
failure are more complex. While the matrix material in the composite may have similarities 
to metal in terms of modes of crack opening, composites have fiber fractures and 
delaminations which affect the redistribution of load. This leads to the conclusion that the 
toughening mechanisms are different in composites and metals. Toughening mechanisms 
are physical phenomena responsible for crack growth resistance. Crack-tip plasticity, for 
example, is a dominant toughening mechanism in metals. The toughening mechanisms in 
composites are due to micro-cracking. Even though the toughening mechanisms are 
fundamentally different for metals and composites, the result for both materials is the 
elimination of the stress singularity and the effect on load redistribution. 

Currendy, much research is being done to characterize the damage tolerance of various 
composite aircraft structures. For example, critical issues surrounding the advanced 
composite fuselage. Figure 2, are the catalysts for many research programs. Consider the 
crown region of a composite fuselage. Hoop and longitudinal stresses in the crown section 
are caused by cabin pressure while additional longitudinal stresses are due to the 
empennage forces on the aircraft during flight. If the fuselage suffers a through penetration 
by some foreign object, the "flaw" or "notch" would reduce the residual strength of the 
structure. Residual strength prediction capability for curved suffened composite structures 
such as a composite fuselage does not yet exist. To develop a progressive damage 
methodology capable of accurate residual strength predictions for a damaged composite 
fuselage, accurate damage modelling of a much smaller scale must be accomplished first. 

As a first step, the behavior of center-crack tension composite laminate panels with 
through penetration notches has been studied. The center-crack tension composite 
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laminate. Figure 3, simulates a small region in the crown portion of the fuselage and 
assumes the through penetration occurs perpendicular to either the hoop stresses or the 
longitudinal stresses. In addition to disregarding the curvature of the fuselage, it narrows 
in on the notch-tip region and assumes the through penetration occurred away from all 
stiffeners. The following literature survey will briefly describe others' contributions in the 
area of damage tolerance and residual strength of composite structures. 

Literature Survey 

Many fracture models are being developed to model damage and predict residual 
strength. Whitney and Nuismer (WN) [3,4] developed the "point-stress" and the "average- 
stress" failure criteria. Both criteria assume fracture occurs when the stress at some 
characteristic distance from the crack tip equals the unnotched strength. Pipes, 

Wetherhold, and Gillespie (PWG) proposed a fracture model to predict the notched 
strength of composite laminates [5]. In this model they claim the characteristic distance 
used in the WN model is not a material constant Finally, the inherent flaw model 
developed by Waddoups, Eisenmann, and Kaminski [6] is an applied classical LEFM 
model which utilizes a characteristic distance and unnotched laminate strength to predict the 
notched laminate strength. 

A thorough study of the fracture toughness and residual strength of various fibrous 
composites was done by Poe [7-11]. Poe found in his investigations of brittle laminated 
composites that linear elastic fracture mechanics (LEFM) could be used to determine the 
fracture toughness of a notched composite panel without loading it to failure. Harris and 
Morris [12-16] conducted a thorough investigation of translaminate fracture in notched 
composite laminates. The influence of stacking sequence is documented as well as the role 
of delamination in thick notched composite laminates. Many observations were made 
concerning fiber fracture, delamination, matrix cracking, and the influence of laminate 
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thickness. Damage was documented using x-ray radiography, and in some cases, 
specimen deply techniques were used. Some successful strength predictions were made 
using LEFM, including Poe's general fracture toughness parameter. Harris [17] also 
conducted an investigation into the use of crack-tip opening displacement with a Dugdale- 
type model to predict notched laminate strength. Similar investigations by Poe et al. [18] 
provides more deplied laminates revealing the ply-by-ply fiber damage. 

Poe et al. [19] continued with damage tolerance studies by considering the crack 
growth resistance of large fuselage panels with through penetrations that represent discrete 
source dam age Crack growth resistance was plotted using the fracture data from various 
tests and residual strength predictions were made using LEFM. Poe found that LEFM 
predictions were too conservative, but crack growth resistance curves (R-curves) 
determined from strengths of unstiffened sheets made reasonably accurate predictions. 
Further discussion of the general fracture toughness parameter is included. Orange [20] 
presents a method by which the crack growth resistance is estimated from residual strength 
data for notched laminates. Information from simple test results can seemingly then be 
used to estimate the failure loads of more complicated structures of the same material and 
thickness. Schwalbe [21] provides an investigation into crack-tip opening displacement 
and crack growth resistance. R-curve methodolgy is explained in detail and driving force 
predictions are made reasonably well. 

The characterization of damaged notched composite laminates has led to numerous 
damage growth models. Aronsson and Backlund [22] used a damage zone model to 
predict strength and load vs. displacement behavior. Damage is represented by a 
Dugdale/Barenblatt cohesive zone where the cohesive stresses decrease linearly with an 
increase in crack opening. Chang et al. [23] used the progressive damage analysis 
developed by Chang [24] to predict damage growth and failure of an open-hole tension 
composite specimen. The failure analysis consisted of property degradation models and 
failure criteria for matrix cracking as well as fiber and fiber-matrix shearing failure. The 
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AUen-Harris model [25-29] is a continuum damage model which utilizes kinematics based 
volume averaged damage variables to represent matrix crack growth and fiber fracture. 

This model has a mode I matrix crack growth law for fatigue as well as monotonic tension. 
As of now, Delamination and Mode II matrix crack growth is modeled empirically. Recent 
works extending from this model can be found in the literature [30-35]. An experimental 
verification of its ability to predict stiffness loss was documented by Coats [33,34] and its 
ability to predict residual strength was documented by Lo et al. [35]. Unfortunately, none 
of the mentioned models provide a means to predict delamination initiation and growth and 
the resulting stress redistribution. There is, however, much work being done in that area. 

Ko et al. [36] predicts delamination initiation load and location by averaging the stress 
over the characteristic length in conjuction with the Hashin-Rotem criterion [37,38]. This 
research was applicable to balanced symmetric laminates containing a hole. Eason and 
Ochoa [39] incorporated a shear deformable theory, to predict out-of-plane shear and 
normal stress, into a finite element formulation for a plate with an open hole loaded in- 
plane. This allowed for the approximation of interlaminar stress, and as a result, the 
prediction of delamination initiation and growth. A technique for calculating strain-energy- 
release rate, G, for delamination around an open hole is developed by O'Brien and Raju 
[40]. The location of delamination around the hole boundary was successfully predicted 
for quasi-isotropic laminates. Lagace and Saeger [41] developed a methodology for 
predicting delamination initiation at holes in composite laminates. They used the 
interlaminar stress state in conjunction with a mechanics of materials failure criterion to 
compute the delamination initiation load. 

All of the previous mentioned works were investigations or modelling with the goal to 
arrive at a prediction methodology that will predict the residual strength of notched 
composite laminates by taking into account all of the failure mechanisms. There are certain 
mechanisms that will redistribute the load at or around a notch. This load redistribution can 
cause a reduction in stress concentration and an increase in strength. Such mechanisms 
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must be accounted for if accurate residual strength predictions are required. For example, 
axial splits at a notch are a means of load redistribution and therefore reduce the stress 
concentrations at the notch. The axial splits, therefore, are known as a toughening 
mechanism. The key is to identify toughening mechanisms and incorporate them into the 
modelling. 

The toughening issue will be discussed thoroughly in Chapter II and Chapter ID. 
Chapter II, Experimental Characterization, will discuss in detail the experimental 
procedures and the experimental verification of crack growth resistance. Since toughening 
mechanisms are responsible for crack growth resistance, experimental evidence and a 
postulate concerning toughening mechanisms is discussed. The Progressive Damage 
Analysis, Chapter III, addresses the issue of modelling the damage. The mechanics for 
matrix cracking and fiber fracture are discussed as well as the mathematical framework for 
the progressive damage model - the Allen/Harris non-linear constitutive model. The 
progressive damage analysis scheme consists of a damage-dependent finite element 
analysis implemented into the NASA Computational Structural Mechanics Testbed 
(COMET). Accounting for the gradual load redistribution effects of toughening is 
discussed in this chapter and analytical results are presented. 

Objectives and Approach 

Since fibers are the major load bearing component in most composites, then predicting 
residual strength of a center-notched composite laminate would require a knowledge of 
fiber fracture as well as matrix crack growth. Therefore, the research herein takes an in- 
depth look at the failure mechanisms involved in translaminate fracture of center-crack 
tension composite panels. Furthermore, fiber bridging in the sense of intact fibers of one 
ply of a laminate bridging the fractured fibers in an adjacent ply is considered as a possible 
explanation for crack growth resistance as stable tearing occurs from the notch-tip. 
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The objective of this study is to develop a progressive damage methodology capable of 
predicting residual strength of composite structures. The approach is as follows: 

• Experimental characterization of failure mechanisms by the following methods: 

(a) x-ray radiography 

(b) ply-level fractography 

(c) R-curve behavior 

• Theoretical postulate resulting from the experimental study: a fiber bridging effect 
is present and is a mechanism of crack growth resistance. 

• Develop fiber fracture failure criteria and implement them into the Allen-Harris 
progressive damage model. 

• Utilize a damage dependent finite element code where the damage is modeled using 
the Allen-Harris model [25-29]. 

(a) mesh refinement study. 

(b) residual strength predictions using a ply discount and an elastic-perfectly 
plastic monotonic damage growth law. 

(c) compare the two damage growth laws to illustrate the bridging effects in the 
elastic-perfectly plastic growth law. 

(d) R-curve predictions. 
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CHAPTER II 

EXPERIMENTAL CHARACTERIZATION 


This Chapter will provide a detailed description of the experimental study. Discussed 
first will be the experimental procedure, i.e. specimen preparation, test setup, and loading. 
The next discussion will be on specimen deply techniques and R-curves to illustrate crack 
growth resistance. Finally, fiber bridging will be introduced as a possible toughening 
mechanism. 

Experimental Procedure 

Material and Specimen Configuration 

Through-penetration of an aircraft fuselage is simulated on a small scale using the 
center-crack tension (cct) specimen. Figure 3. Five configuration groups were tested; 1" 
unnotched coupons for collecting unnotched material properties, 4" wide and 12" wide 
specimens where the 4" wide specimens had 1/2" notches (2ao=l/2", w=4", L=18") and 1" 
notches (2ao=l", w=4", L=18"), the 12" wide specimens had a 3" notch (2ao=3", w=12", 
L=34"), and the 36" wide panels had 9" notches (2ao=9\ w=36", L=90"). Each 
configuration group consisted of the materials AS4/8553-40, AS4/938, and AS4/3501-6. 
The AS4/3501-6 was manufactured using a tape pre-preg while the others were made from 
the tow-placement technique. The AS4/938 specimens had two different tow spacings. All 
of the test materials with their test identification name, material identification, dimensions, 
and manufacturing techniques are summarized in Table 1. The layups used were 
[+ 45/0/90/+ 30/d]s for the center-crack tension specimens and its transverse 


8 



[+45/45/90/0/^60/ 90] s was used for some of the unnotched coupons as well. The 
specimens were fabricated by Boeing and shipped to NASA Langley Research Center for 
testing. The laminate stacking sequences are for fuselage structures designed by Boeing. 

Specimen Preparation 

Each specimen was strain gaged according to Figures 4-7c. The 1" unnotched coupons 
and the 4"wide panels were monotonically loaded to failure in a 50 kip servo-hydraulic 
testing machine. The 12" and 36" wide panels had anti-buckling guide plates attached just 
above and below the notch. They were loaded to failure in a 100 kip and 500 kip servo- 
hydraulic testing machine, respectively. The strain gages were all wired to a Vishay 
Measurements Group System 4000 data acquisition unit. The system 4000 also collected 
the applied load from the load cell and the center-crack opening displacement from the ring 
gage secured in the center of the notch. 

Loading and Data Collection 

As the panels were being loaded, the discrete source damage (fiber fracture, 
delamination, and matrix cracking local to the notch tip) was frequently audible. 
Periodically, as damage progressed with increasing load, a zinc-iodide dye penetrant was 
applied to the notch and edge of the specimen. X-ray radiographs were taken of the right 
and left notch-tip regions. The damage absorbed the zinc-iodide dye penetrant and the 
damage is accurately represented in the x-ray radiograph as a blackened or shaded region. 
Figure 8. Lamina material properties were obtained from the literature and are provided in 
Table 2. The laminate properties including failure loads are given in Tables 3 thru 8. The 
accuracy of the lamina material properties is somewhat questionable since the lamina 
material properties were obtained from various references and it was nearly impossible to 
find a complete and consistent set of data. Lamina properties for AS4/3501-6 were taken 
from an ASTM STP [42] and the other properties were taken from Boeing test data [43]. 
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Some of the panels were not loaded to catastrophic failure. Instead, x-rays were taken 
periodically up to a percentage of the ultimate failure load. The specimen was then taken 
out of the grips and the area surrounding the notch tip was deplied. Implications and 
conclusions about crack growth resistance from the x-ray radiographs and specimen deply 
techniques are discussed in the following paragraphs. 

Verification of Crack Growth Resistance 

There are a number of ways to verify crack growth resistance. To do this, an 
understanding of what crack growth resistance means is necessary. Crack growth 
resistance is simply a resistance to crack propagation. A damage tolerant material has a 
high resistance to crack growth and the physics of this resistance can be illustrated in x-ray 
radiographs of the notch-tip damage, fractographs of deplied specimens, and R-curves. 

X-Rav Radiography 

Typical damage tolerance as might be seen in a composite fuselage can be illustrated in 
an x-ray radiograph of a center-crack tension panel. Figure 8. Notice in this figure the 
amount of damage accumulation, or the damage tolerated. The left and right notch tip 
damage is shown. What looks like a "tear" in the panel is fiber fracture, and the shading 
surrounding it is local delamination. The lines extending from the notch tip and through the 
fractured region are matrix cracks in the off-axis plies. The key issue here is the amount of 
crack growth before catastrophic failure. This x-ray radiograph is evidence of crack 
growth resistance in the composite fuselage because the panel is still sustaining the very 
load that created such damage. Notice that in Figure 8 the damage seen is only at 89.6% of 
failure when this x-ray radiograph was taken. This is an indication of a damage tolerant 
structure. 
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A series of x-ray radiographs from various size panels are provided in Figures 9-1 lb to 
illustrate crack growth resistance. It was deemed unnecessary to provide the x-ray 
radiographs for every specimen since damage initiation and growth occurred similarly for 
all of the specimens anyway. Notice the loads at which these x-ray radiographs were 
taken. Crack growth resistance is what keeps the fiber fracture from traversing the entire 
width of the laminate. The mechanism of such resistance will be discussed later in this 
chapter in Isolation of Toughening Mechanisms. The point here is that because the x-ray 
radiograph shows such extreme damage at only a percentage of ultimate failure, crack 
growth resistance exists. 

Fractograohy of Pl v Level Damage 

A few specimens were chosen to be loaded only up to a percentage of ultimate failure. 
Then the notch-tip damage region was isolated and pyrolized in an oven at 850°F until the 
neat resin had burned away (about four hours). The individual plies were separated and 
examined using a JEOL JSM-6400 scanning electron microscope. Fiber fracture is clearly 
visible in these plies. Figures 1 2a-c. The zinc iodide stain produced from the x-ray 
procedure is a reliable indication of local delamination. A schematic of the fiber fracture 
and local delamination is provided in Figure 1 3 to aid the eye in locating and quantifying 
the ply damage. The dark patches drawn on the schematics at the notch-tips represent local 
delamination and fiber fracture is represented by the "free-hand” drawn lines. Dimensions 
are given for most of the delaminations and the fiber fracture is dimensioned as da. 
Schematic representations of more fractographs are provided in Figures A-l to A-3 of 
Appendix A. The schematics illustrate the evidence of crack growth resistance as well as 
the fractographs. It is not necessary to show all of the fractographs and schematics because 
they are repetitive illustrations of the various patterns and magnitudes of fiber fracture in 
each ply. Therefore, the conclusion that crack growth resistance is clearly evident is the 
same for all of the fractographic examinations. 
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R-Curve Behavior 

During the monotonic loading of these panels, the applied load and the crack opening 
displacement (cod) at the center of the notch was recorded and used to produce load/cod 
plots. Figures 14-16. Discontinuities, or jumps, exist at various places along the load/cod 
plot where the fiber fracture was audible during loading. At these discontinuities the 
specimen was unloaded to take an x-ray and are labeled A, B, C, etc., on the load/cod plot. 
The corresponding x-ray radiographs are given in the plots to illustrate the amount of 
damage at each discontinuity. A closed form elasticity solution [44] for determing the 
characteristic half crack length, a, for a quasi-isotropic material under plane stress is 


COD = 


4*S»a 


0 ) 


where Ex is the longitudinal modulus, S is the applied stress, and a is the characteristic half 
crack length. Given that Aa=a-ao, and after some algebraic manipulation. 


Aa = a 0 


W ■ ' ) 


(2) 


is the elasticity solution for effective crack growth and is shown in the figures as well. The 
initial half crack length is ao, and the subscript i indicates the initial load/cod slope up to the 
point of separation between the load/cod curve and its initial slope. The load/cod plots 
illustrate the use of this closed form solution. 

A plot of fracture toughness as a function of effective crack growth is called a crack 
growth resistance curve (R-curve). The fracture toughness, in terms of the general fracture 
toughness ratio, of large notched composite laminates is given in Figure 17 [19]. The 
general fracture toughness ratio, Qc/£tuf> was developed by Poe [7-11] using linear elastic 
fracture mechanics (LEFM) where Q c is the general fracture toughness parameter and 8tuf 
is the tensile failing strain of the fibers. Qc is independent of laminate orientation and was 
derived on the basis of fiber failure in the principal load-carrying laminae. Q c is 


12 



proportional to the critical value of the mode I stress intensity factor (S.I.F.) and the 
constant of proportionality depends only on the elastic constants of the laminate. Poe 
showed that the ratio Qc/£tuf was a constant for all brittle epoxy composite laminates 
regardless of layup. Therefore, a single value of Qc/£tuf could be used to predict the 
fracture toughness of these fibrous composite laminates from only the elastic constants and 
8tuf. Experimental data [7-10] indicated Qc/£tuf is reasonably constant, 0.3 Vin , except 
for instances where extensive delamination or 0° ply splitting occurred. 

A couple of the laminates in Figure 17 exhibited significantly higher fracture toughness 
than Poe's prediction. Furthermore, the crack growth resistance evident in Figure 18 [19] 
portrays a much higher fracture toughness than Poe's previous investigations. This is due 
to toughening mechanisms not accounted for in the constant general fracture toughness 
parameter. The various toughening mechanisms affect crack growth resistance and thus 
affect the shape of the R-curve and the value of the fracture toughness. 

The effective crack growth calculated from the load/cod plots mentioned above and the 
corresponding applied stress is used in Poe's general fracture toughness solution [7-11] 

S T , . „ Tc(a+Aa)"|l/2 K 

Q c = 1 -[K(a + Aa)sec — ^ — J (Q ( 3 ) 

to generate crack growth resistance curves (R-curves) in Figures 19-21. An x-ray 
radiograph for one of the specimens in each graph is supplied to illustrate the extent of the 
notch tip damage. The dashed line noted as previous work is the constant Qc/£tuf 
determined in Poe's previous investigations [7-11] and is placed in the figures to illustrate 
that the fracture toughness of some notched composite laminates is not accurately predicted 
by the constant ratio. 

The R-curves are experimental verification of crack growth resistance because they 
illustrate the continuing load carrying capability with increasing discrete source damage. 
From the three experimental verification techniques, x-ray radiography, scanning electron 
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microscope fractography, and R-curve behavior, we can state that these center-crack 
tension panels were damage tolerant since much of the damage was detectable long before 
catastrophic failure. The crack growth resistance was due to a dominant toughening 
mechanism and these toughening mechanisms need to be accounted for in a progressive 
damage analysis if accurate residual strength predictions are to be obtained. 

Isolation of Toughening Mechanisms 

Ample evidence of crack growth resistance was presented in the previous sections. 
Materials that are damage tolerant resist crack growth because of one or more existing 
toughening mechanisms. It has been proposed and accepted by many researchers that fiber 
bridging, intact fibers bridging the wake of a matrix crack, is a dominant toughening 
mechanism in many materials. Fibers bridging matrix cracks is not a likely or realistic 
toughening mechanism for the center-crack tension panels investigated in this study. 
However, we may postulate that fiber bridging in the sense of intact fibers of one ply 
bridging the fractured fibers of another ply is a dominant toughening mechanism by which 
load is redistributed. It is obvious from the x-ray radiographs in Figures 9-11 and in the 
ply fractographs in Figures 12a-c that this type of fiber bridging is an existent physical 
phenomenon in the center-crack tension composite panels. Throughout the remainder of 
this study, this type of bridging will be referred to as ply bridging. (Appendix B provides 
a discussion on fiber bridging and ply bridging, including a literature survey of fiber 
bridging). The next step involves applying the load redistribution effect of the ply bridging 
in a progressive damage analysis so that it will be possible to predict residual strengths of 
the center-crack tension panels within an acceptable level of accuracy. 
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CHAPTER III 

PROGRESSIVE DAMAGE ANALYSIS 

Fiber fracture criteria were developed and implemented into an existing progressive 
damage model framework, the Allen-Harris model [25-29], The progressive damage 
model was implemented into a multi-purpose finite element code [35] and a residual 
strength prediction capability was developed. This progressive damage methodology is 
damage dependent and can therefore model the damage development at and around the 
notch-tip. It is also independent of laminate stacking sequence, and the finite element 
analysis makes it possible to analyze any geometrical configuration. 

A Progressive Damage Model 

The Allen-Harris Model 

The damage model of Allen and Harris [ 25-29] was originally developed to model the 
behavior of microcrack damage in brittle epoxy systems and has recently been extended to 
toughened polymer systems. The model predicts the growth of intraply matrix cracks for 
monotonic tensile loadings and for tension-tension fatigue, the associated ply level damage- 
dependent stress and strain states, and the residual strength of laminates with geometric 
discontinuities. The model also accounts for the effects of delaminations but uses an 
empirical relationship that requires the user to supply an estimate of the delamination area. 
The empirical relationship must be used because the model currently does not calculate free 
edge interlaminar stresses. (The mathematical formulation of the model may be found in 
the literature [29] and will not be reproduced herein.) The model uses internal state 
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variables (ISV) to represent the local deformation effects of the various modes of damage. 
Loading history dependence is modelled by ISV damage growth laws. The progression of 
damage is predicted by an iterative and incremental procedure outlined in the flowchart 
shown in Figure 22. This entire progressive failure analysis scheme has been implemented 
into the finite element formulation in the NASA Computational Mechanics Testbed 
(COMET) [35] computer code. The first block of Figure 22 is a description of the 
information needed as model input. A FORTRAN code consisting of the damage 
dependent constitutive model and a damage growth law for matrix cracking was 
incorporated into a classical lamination theory analysis to produce effective lamina and 
laminate properties for unnotched laminates. The program is called FLAMSTR (Fatigue 
LAMinate STRess) [32] and makes up the first constitutive module. The fourth block is a 
damage dependent finite element analysis code [35] from which the second constitutive 
module performs a ply level elemental stress analysis and simulates damage growth via 
damage growth laws for each element. The damage growth calculations, block six, are 
used to update the damage state, block seven, for the notched laminates. Note that for 
unnotched laminates, only the first constitutive module is needed to update the damage 
state. 

The material property descriptions required for the model include standard ply stiffness 
and strength data determined in the usual manner. In addition, the tension-tension fatigue 
matrix crack growth law must be determined from test data obtained from the [0/902/0] s 
laminate. Under tension-tension fatigue, matrix cracks accumulate in the 90 degree layers 
and, therefore, the effects of mode I mauix crack growth is isolated. The mode II matrix 
crack growth law can be obtained from fatigue tests of the [45/-45] s laminate which isolates 
the 45 degree plies in pure shear. (The mode II growth law is not currently implemented 
into the finite element code.) A procedure [33] has been developed for determining the ISV 
(damage parameters) from the test data obtained from these two laminates. 
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Quasi-Static Loads and Damage Modelling 

Recall that the model uses internal state variables to represent the local deformation 

effects of the various modes of damage. For instance, the ISV representing mode I matrix 

cracking is of* . When the material is subjected to quasi-static (monotonic) loads, the rate 
L 22 

of change of the internal state variable, ^ orm 

, M _ fPd(£22-e22crit) if £22 > £22crit; ^ 

L 22 lO if £22 < £22crit 

where £ 22 crit is the critical tensile failure strain and P is a factor that describes the load 
carrying capability of the material after the critical tensile strain has been reached. A similar 

relationship is used to describe the tensile failure of the reinforcing fibers. The internal 

state variable for this mode of damage is a, M and its rate of change is of the form 

Ml 

i M _ fyd(£| l~£l lerit) if ^11 ^ ^llcrit* 

^11 "lo if e 1 1 < Client 

where £ii crit is the tensile fiber fracture strain and y is a factor describing the residual load 
carrying capability of the material after fiber fracture has occurred. The numerical details of 
0° ply fiber fracture as they appear in the finite element code are as follows. The 
longitudinal 0° ply stress (ignoring thermal strains) is written as 


c 1 =Qi,[e 1 < w ]+Q, 2 [e 2 -a* w ] (6) 


where ai is the longitudinal 0° ply stress, Qn and Q 12 are the ply level reduced moduli, Ej 
and £2 are the 0° ply longitudinal and transverse strains, respectively, and a" ew and a^ ew 

represent the updated 0° ply local deformation effects of fiber fracture and matrix cracking, 
respectively. The change in the ISV representing the effects due to tensile fiber fracture is 


doti = 



a 0 /* 1 ) + Q12 (e 2 - a° 2 ,d ) ' Y s x 
Q11 


(7) 
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where is the lamina longitudinal critical strength. The ISVs are updated by 


a*" = a° Id + don (8) 

and then equation (8) is substituted into equation (6). The result of this substitution is 

0i = Qi i [e i - | cij ld + d ot i | J + Qi2^£2 - (9) 


which is further modified by substituting equation (7) into equation (9) to obtain 


Gi = Q 


n 


o ~old 

El'Ot j 


|ei-a 


old . Qx2 

Qn 


(£2-a 2 ) - 


QnJ. 


+ Ql2 (£2-012) 


( 10 ) 


and 

0i=7Sy (11) 

is the result of obvious cancellations. Notice that if y=0, the monotonic failure criterion 
results in a ply discount type of behavior. If 7^=1, elastic-perfectly plastic behavior is 
obtained. A computer algorithm has been written for this computational scheme and 
implemented into a finite element analysis code. 

Finite Element Analysis 

Model Configuration and Mesh Refinement 

A previous mesh refinement study showed that the analytical solutions for residual 
strength converges very well for open-hole laminates. (The open-hole mesh refinement 
study is documented in Appendix C). In the analysis of the center-crack tension panels, an 
initial mesh refinement study revealed that a very fine mesh in the notch region severely 
under-predicted the failing load, and a course mesh resulted in a failure load much higher 
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than the experimental values. The fine mesh analysis revealed out-of-plane displacements 
and in-plane rotations were occurring at and around the notch-tip on the order of 10' 1 . The 
course mesh results showed no such in-plane rotations or out-of-plane displacements more 
than on the order of 10' 14 . It was decided that these rotations and displacements must be 
an indication of localized buckling. 

The reliability in the experimental and analytical strain field correlation was previously 
demonstrated by the stiffness loss predictions for the open-hole laminates [34]. Therefore, 
a comparison of the model predicted buckling effect to the experimental strain gage data 
was observed to confirm the existence of localized buckling. Typical notch-tip stress/strain 
behavior for the center-crack tension panels. Figure 23, is an obvious illustration of 
localized buckling. This dcmonsuates that the experimental stress/strain behavior correlates 
well with the analytical displacements and rotations to conclude that localized buckling is 
indeed occurring at and around the notch-tip of the center-crack tension panels. This is a 
reasonable conclusion since localized buckling in center-crack panels under tensile loadings 
is a well known phenomenon. Sawicki et al. [45] documented this phenomenon in their 
photoelastic investigations of center-crack tension panels. 

The mesh refinement study also addressed the issue of choosing an optimum mesh that 
provides the best results. The model averages the kinematic effect of damage over the 
entire element. This averaging procedure results in a length scale which is an important 
consideration in mesh generation. Since the progressive damage model represents damage 
with volume averaged quantities that are averaged over the entire element, an element too 
large may not represent the effects of the stress singularity and results in an over-prediction 
of strength. If the element area relative to the notch size and the material constituents is 
small, that element size approaches the micro scale. An element too small may cause the 
averaging process to exaggerate the effects of the stress singularity and result in an under- 
prediction of strength. This is similar to the Whitney/Nuismer point stress or average 
stress criteria [3,4], Based on experimental data, there seems to be a characteristic distance 
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which determines the proper element size. This phenomenon is also not unique to 
composites. There has already been much discussion on this topic relating to crack-tip 
plasticity in metals and dates back at least as far as the investigations by Chan et. al. [46]. 

The finite element model configuration is a quarter panel mesh of a center-crack tension 
panel with a 0.02" notch. Figures 24-27. These meshes have an element size and spacing 
at the notch-tip which resulted from the combined mesh refinement and localized buckling 
study. The meshes were constructed so as to allow the existence of localized buckling and 
obtain reasonable analytical solutions. All of the nodes on the y-axis midplane from the 
notch tip to the panel's edge are constrained in the x-direction, the nodes on the x-axis 
midplane are constrained in the y-direction, and all the nodes are constrained in out of plane 
rotation. The loading is applied in the x-direction. 

Analytical Predictions 

Residual strength predictions have been made for center-crack tension laminates loaded 
in monotonic tension. R-curves for the center-crack panels were predicted and all 
predictions were compared with experimental data. These results will be discussed in the 
next two sections. 

Residual Strength Predictions 

The progressive damage model computed residual strengths using two different failure 
criteria (monotonic damage growth laws). The first law, a ply discount criterion, was 
achieved by setting the monotonic growth law parameter, y, to zero (T=0). Recall from 
equation (11) that when fiber fracture occurred, if y=0, the load carrying capability of a ply 
within an element would be eliminated. This criterion does not account for any toughening 
mechanisms and is therefore extremely conservative. The second law, achieved with y^l, 
is an elastic-perfectly plastic criterion. The longitudinal ply stress can only be as high as 
the ply critical strength of the ply at or after fiber fracture in this case. This criterion 
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accounts for load redistribution and is more likely to represent some of the toughening 
mechanisms seen in the experiments such as ply bridging. 

The 0° ply longitudinal stresses were plotted along the transverse distance from the 
notch-tip for the AS4/938 panel with the 3" notch. The plots for both failure criteria are 
plotted together prior to any damage in Figure 28 to demonstrate that for no damage, the 
stress calculations are consistent and the failure criteria has no effect prior to damage. Ply 
discount is illustrated in Figure 29 by the portion of the plot that shows a ply stress value of 
zero psi. Fracture in the 0° plies has occurred in the elements at a transverse distance away 
from the notch-tip at almost 0.9". Along this distance the load carrying capability of the 
fractured plies has reduced to zero. The load carrying capability plateaus at the ply critical 
strength for the y=l case in Figure 30. There are load cycles when the notch-tip stresses 
reach negative values for both cases y=0 and y= 1 . This is due to the extreme displacements 
and rotations occurring at the notch-tip. The change in the internal state variables is not as 
drastic as the strains in this region and as fracture progresses away from the notch, the 
notch-tip strains decrease, the internal state variables do not decrease, and the resulting 
stress is negative. This is a numerical artifact in the code that will be corrected in the near 
future. This does not affect the residual strength results. 

The residual strength predictions for y=0 and y=l are illustrated in Figures 31-33 where 
the model predictions are compared to the experimental averages. The error bars represent 
the experimental minimums and maximums when available. The experimental values were 
given in Tables 4-8 and the values represented in Figures 31-33 are given in Tables 9-11. 

The predictions were reasonable considering the model depends on accurate material 
properties. The lamina material properties used in this investigation were not 
experimentally measured from a sample of the material tested. The literature was searched 
for a consistent set of material properties for AS4/3501-6 and the only complete set found 
was in a journal article [42], The other properties came from data documentation by 
Boeing [43]. 
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A lamina material properties sensitivity study revealed that a change in only one 
property could have a noticeable affect on the analytical solution. In Tables 12 and 13, 
three analytical studies were performed. The first analytical column provides the actual 
chosen properties for this study and the resulting solution is listed at the bottom of the 
table. The second and third analytical columns each have a property in italic and bold faced 
font These values were the only change for that analysis and the resulting solutions are 
given at the bottom of the tables. The experimental column has the experimental minimum 
and maximum residual strength. Notice that for AS4/8553-40, a 1.5% increase in En 
results in a 1.8% increase in strength. Furthermore, an 1 1.5% increase in critical strain 
results in a 4.7% increase in residual strength. For AS4/3501-6, a 10% decrease in E22 
had no effect as expected, and a 10% decrease in Gi2 resulted in a 4.7% decrease in 
residual strength. No attempt was made herein to match the experimental results by 
selecting material properties that "optimized" the predictions. 

The residual strength predictions are not only sensitive to the lamina material properties 
but to the fiber fracture failure criteria as well. The decision to choose y=l (elastic-perfectly 
plastic) and y=0 (ply discount) was not a random thought or an attempt to find a criterion 
that would match experimental results. The school of thought here was to have two 
extremes that would allow a study of the toughening mechanism caused by the ply bridging 
without any other phenomenological influences that would fit the analytical solutions to the 
experimental data. Both of the failure criteria allow load redistribution through ply 
bridging. Unlike the elastic-perfectly plastic criterion, the ply discount method does not 
allow a fractured ply in a given element to carry any load. The load redistribution is sudden 
and the adjacent plies fail very quickly after the first ply failure. This is why the ply 
discount method is said not to have any toughening mechanisms even though load 
redistribution is occurring. 


22 



If the ply discount method was able to predict the residual strengths of the center-crack 
tension panels within +10%, there would be no need for this progressive damage model or 
a fiber fracture failure criterion that accounts for the toughening effects observed in the 
experiments. The model, utilizing the elastic-perfectly plastic fiber fracture failure criterion, 
made predictions that fell within +10% of the experimental averages in most cases. Most 
of these predictions fell within or very near the experimental data scatter. Considering that 
the material properties may vary +10% depending on where and how they are obtained, the 
residual strength predictions are about as good as they can get with the exceptions of 
perhaps the 1" notch of the AS4/8553-40 and the 9" notch of the AS4/3501-6. It was 
believed that the failure load of the 9" notch AS4/3501-6 was much lower than it should 
have been given the lamina material properties of this material. The lamina material 
properties for the AS4/8553-40 was found in a NASA Contractor Report [43]. The 
properties was for a material with a fiber volume fraction of about 53% and the actual 
material tested in this study had a fiber volume fraction of 58%. The rule of mixtures was 
applied to obtain the lamina material properties used in the model, and there is obviously 
compounding errors due to that step as well. 

These results show that the model has reached a level of maturity where it can be used 
to model fiber and matrix damage progression and predict the residual strength of notched 
composite laminates. This study was key in developing the model to this level of maturity. 
The predicted R-curves will demonstrate this further and conclude the analytical results of 
this investigation. 

R-Cur\>e Behavior 

The R-curves are plotted in Figures 34-45 using the elastic-perfectly plastic monotonic 
damage growth law. These figures show model generated R-curves compared to the 
experimental R-curves. The effective crack growth, Aa, is calculated using the closed form 
solution in equation (2). R-curves were generated using actual Aa measurements from the 
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x-ray radiographs in Appendix D. However, due to the subjectivity of measuring the 
damage from an x-ray radiograph, the closed form solution approach in Figures 34-45 will 
be the center of discussion in this section. 

Notice in these figures that the general shape of the model generated R-curve is similar 
to the experimental R-curves except for the initial data points. These initial points of 
fracture are critical because they determine the slope and curvature of the initial portion of 
the curve. The R-curves reveal that the predicted initial fiber fracture occurs at a higher 
fracture toughness than the experimental data. However, as loading continued, an under- 
prediction of residual strength was manifested in a lower fracture toughness in the model 
generated R-curves. Likewise, the model generated R-curve revealed a higher fracture 
toughness than the experimental fracture toughness where the residual strengths were over- 
predicted. There was one exception in Figure 46 where it would seem only panel 
G2TAPEA failed at a higher load than the model prediction, when in fact panel F1TAPEA 
had a higher residual strength as well. The problem with modelling the initial fiber fracture 
at a higher fracture toughness is evident here and is a characteristic of the monotonic 
damage growth law, lamina material properties, and the critical failing strains used in the 
model. 
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CHAPTER IV 

CONCLUDING REMARKS 


An experimental investigation of translaminate fracture was conducted and a residual 
strength prediction capability was developed using a progressive damage methodology. 
An experimental characterization of several composite materials systems revealed an R- 
curve type of behavior. Fractographic examinations led to the postulate that this crack 
growth resistance could be due to ply bridging, defined in this study as fractured fibers of 
one ply bridged by intact fibers of an adjacent ply. The Allen-Harris model was used in a 
finite element code to model the matrix cracking and fiber fracture that results from the 
notch-tip stresses in center-crack tension composites. Two fiber failure criteria were used 
to model the progression of fiber fracture. The first criterion is essentially the classical 
ply discount method because as fiber fracture occurred in a particular ply for any given 
element, that ply in the given element would no longer have any load carrying capability. 
Therefore, this criterion did not allow for any toughening effects. The other fiber failure 
criterion is an elastic-perfectly plastic fiber failure growth law. As a ply fractures in an 
element, this fiber failure criterion allows for a more gradual load redistribution and the 
load carrying capability of that ply in the element is constrained to the lamina 
longitudinal failure strain. These criteria were chosen as two extremes to compare the 
results and obtain an understanding of the ply bridging effect in the elastic-perfectly 
plastic criterion. Perhaps a criterion somewhere in between that is more representative of 
strain softening [47] would provide more accurate results. 

Residual strength was predicted using both of the fiber failure criteria. The elastic- 
perfectly plastic criterion resulted in predictions within + 10% of the experimental 
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averages in most cases. Furthermore, because the elastic -perfectly plastic criterion is 
more representative of the toughening mechanisms observed in the experiments, the ply 
bridging concept was shown to be valid for predicting residual strength for tensile 
dominated failure loads. Predictions of the R-curve type of behavior were made as well. 
These predictions were relatively consistent with the residual strength predictions in that 
under-predicted residual strengths also resulted in under-predicted R-curves. 

This investigation was only a small step in the efforts to develop a progressive 
damage methodology to predict residual strengths of composite aircraft structures. For 
this particular model, not all of the mechanisms of damage have been included. 
Delamination initiation and growth still needs to be incorporated into this model as well 
as compression damage mechanisms. All of these mechanisms contribute to the failure 
process as well as load redistribution. It is difficult to say, for instance, whether or not 
the residual strength predictions will decrease if the mechanism of delamination is 
introduced into the modelling. One thought is that delamination would weaken the 
laminate causing a reduction in the residual strength. However, local delamination at or 
around a notch could relieve the high stress concentration and thus increase the residual 
strength. 
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Table 1 Summary of Experimental Test Matrix 
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[+45/0/90/+30/0] s d 0.025" - 0.015" gaps with stacked hands 

[ +45/90/0/^*60/90] s e 0.015" + 0.015" gaps with offset bands 

0.015" + 0.015" gaps with stacked bands 
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Table 3 Data From 1" Wide Unnotched [4-45/0/90/+30/0] s Laminates 
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AS4/8553-40 tow placed (0.015" - 0.015" gaps with stacked bands) 

AS4/938 low placed (0.025" * 0.015" gaps with stacked bands) 

AS4/938 tow placed (0.015" - 0.015" gaps with offset bands) 
AS4/3501-6 tape prepreg 




Table 4 Data From 1" Wide Unnotched [-45/90/0/- 60/90]s Laminates 
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AS4/938 tow placed (0.025" * 0.015" gaps with stacked bands) 

AS4/938 low placed (0.015" + 0.015" gaps with offset bands) 
AS4/3501-6 tape prepreg 










Table 5 Data From 4" Wide Notched [445/0/90/430/0] s Panels (2a 0 =l/2") 
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Table 6 Data From 4" Wide Notched [T45/0/90/*30/0] s Panels (2a 0 =l M ) 
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Table 7 Data From 12" Wide Notched [T-45/0/90/+30/0] s Panels (2a 0 =3") 
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Table 8 Data From 36" Wide Notched [F45/0/90/4-30/0]s Panels (2ao=9 M ) 
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Table 13 - AS4/3501-6 Material Property Sensitivity Stud 
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Figure 1 - In-Flight Structural Failure. 
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Figure 2 - Fuselage Critical Technology Issues. 
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Figure 5a - Strain Gage Schematic for a 4" Wide Notched Panel. 



B( ) - Use EA-06-070LC-350 single gages on front and 
back. Write the gage number in parentheses 
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B(5,6) 

Figure 5b - Strain Gage Schematic at the Notch for the 4" Wide Panel. 


A( ) - Use CEA-06-187UW-350 single gages on front and back. 
Write the gage number in parentheses next to the gage. 



Figure 6a - Strain Gage Schematic for a 12" Wide Notched Panel. 
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Figure 6b - Strain Gage Schematic at the Notch for the 12" Wide Panel. 




Figure 7a - Strain Gage Schematic for a 36" Wide Notched Panel. 
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Figure 7b - Strain Gage Schematic at the Notch for the 36" Wide Pane 



E( ) - Use C-891 1 1 3-C single gages on front and 

back. Write the gage number in parentheses 
next to the gage. 
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Figure 7c - Strain Gage Schematic at the Notch for the 36" Wide 



AS4/938 (;45/0/90/' + 30/0) 
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Figure 8 - Typical Crack-Tip Damage Due to Hoop Stresses in an AS4/938 
Crown Quadrant of a Fuselage Laminate. 
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Figure 9 - X-Ray Radiographs of the Right Notch Tip Damage for the AS4/8553-40 Panels. 
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Figure 10a - X-Ray Radiographs of the Right Notch Tip Damage for the AS4/938 Panels. 
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Figure 10b - X-Ray Radiographs of the Left Notch Tip Damage for the AS4/938 Panels. 
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Figure 1 la - X-Ray Radiographs of the Right Notch Tip Damage for the AS4/3501-6 Panels. 
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Figure 1 lb - X-Ray Radiographs of the Left Notch Tip Damage for the AS4/938 Panels. 
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Figure 12a - Ply Level Fractography of [-45/45/0/90/-30/30/0] s Specimen F5 AK5B Right Notch Tip at 84% of S u jj- 
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Figure 12b - Ply Level Fractography of [-45/45/0/90/-30/30/0] s Specimen F5 AK5B Right Notch Tip at 84% of S u j t - 
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Figure 12c - Ply Level Fractography of [-45/45/0/90/-30/30/0] s Specimen F5 AK.5B Right Notch Tip at 84% of S u jf 
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90 degree ply 




- Fiber Fracture and Delaminationjllustration for Specimen 
F5AK5B, [-45/45/0/90/- 30/30/0] s . Right Notch Tip. 


66 




67 




68 





c 

O 

O 

O 


c 

a> 

E 

a> 

a 

iS 

Q. 

</) 

5 


U) 

c 

’E 

a> 

a. 

O 

ic 

o 

as 

k. 

O 


t? S' 



< 

pa 

Cu 

< 

H 

<N 

o 


0 ) 

c 

c 3 

Oh 

vo 


o 

un 

m 

< 

Uh 

O 


c 

o 

6 

aj 

o 

OH 

a 

c/: 

5 


U) 

c 

c 

a-> 

D- 

o 


ic 

a 

03 

Uh 

u 


vO 


<u 

u- 

3 

W) 

E 


69 





CO CD CM o 

• ■ • • 

o o o o 


O g ? 

O w- C 


70 


Figure 17 - Fracture Toughness of Large Notched Composite Laminates. 
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Effective crack growth, Aa, in. 

Figure 18 - R-Curves for Large Notched Composite Laminates. 
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Crack Growth, Aa (in.) 

Figure 19 - Crack Growth Resistance in AS4/8553-40 Panels. 




AS4/938 
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Figure 20 - Crack Growth Resistance in AS4/938 Panels. 




AS4/3501 -6 [-45/45/0/90/-30/30/0] 
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Crack Growth, Aa (in.) 

Figure 21 - Crack Growth Resistance in AS4/3501-6 Panels. 





Figure 22 - Progressive Failure Analysis Scheme. 
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Figure 24 - Quarter Panel Mesh of the Center-Crack Tension Panel. 

0.02" Notch, W=4", 2a 0 =l/2", L=12", 570 Nodes, 520 
Elements. 
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Figure 26 - Quai 
0.02 



1 00000 
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Distance From Notch Tip, in. 

Figure 28 - 0 Degree Ply Longitudinal Stress vs. Distance from the Notch-Tip. 




82 


Distance From Notch Tip, in. 

Figure 29-0 Degree Ply Longitudinal Stress vs. Distance from the Notch-Tip. 
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Distance From Notch Tip, in. 

Figure 30 - 0 Degree Ply Longitudinal Stress vs. Distance from the Notch-Tip. 
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Figure 31 - Residual Strengths of AS4/938 Center-Crack Tension Panels. 
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Figure 32 - Residual Strengths of AS4/8553-40 Center-Crack Tension Panels. 
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Figure 33 - Residual Strength Predictions of AS4/3501-6 Center-Crack Tension Panels. 
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0.05 0.1 0.15 

Crack Growth, Aa (in.) 

Figure 34 - Experimental and Model Generated R-Curves for AS4/938 
Center-Crack Tension Panels, 2a =1/2”. 
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Center-Crack Tension Panels, 2a = 




AS4/938 
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Center-Crack Tension Panels, 2a 




AS4/938 
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Center-Crack Tension Panels, 2a = 




AS4/8553-40 
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Crack Growth, Aa (in.) 

Figure 38 - Experimental and Model Generated R-Curves for AS4/8553-40 
Center-Crack Tension Panels, 2a =1/2". 
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Center-Crack Tension Panels, 2a = 




AS4/8553-40 
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Center-Crack Tension Panels, 2a = 
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Crack Growth, Aa (in.) 

Figure 42 - Experimental and Model Generated R-Curves for AS4/3501-6 
Center-Crack Tension Panels, 2a =1/2". 
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Center-Crack Tension Panels, 2a = 




AS4/3501-6 
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Center-Crack Tension Panels, 2a = 




AS4/3501 -6 



98 


Center-Crack Tension Panels, 2a = 




APPENDIX A 

PLY FRACTOGRAPHY SCHEMATICS 

Ply fractography for specimen F5AK5B was presented in Figures 12 and 13. The 
schematics for the rest of the ply fractography specimens will be presented in this 
appendix. The actual fractographs will not be shown here since the schematics are easier 
to read. The schematics in Figures A-l to A-3 show representations of local delamination 
as dark patches and fiber fracture as "free-hand" drawn lines. Most of the delaminations 
are dimensioned and the fiber fractures are measured and labeled as da. For example, in 
Figure A-l the first schematic shows the local delamination at the -45/45 ply interface at 
the right notch-tip. The ply used for this was the 45 degree ply. The delamination 
measured 0.31" away from the notch-tip at it's farthest point away. The next ply, a 0 
degree ply had delamination and fiber fracture. The fiber fracture is about 0.20" long 
measured from the notch-tip. All of the dimensions were measured using a machinists 
scale with 1/100 of an inch divisions. All of the figures in this appendix illustrate the 
variation in damage occurring from one ply to another. They demonstrate that as 
translaminate fracture occurrs, the extent of fiber fracture in one ply is not necessarily the 
same as an adjacent ply. This allows for various avenues of load redistribution and in 
effect, existence of toughening mechanisms. 
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45 degree ply 


0 degree ply 


90 degree ply 



90 degree ply 





Figure A-l - AS4/3501-6 Fiber Fracture and Delamination Illustration, 

Specimen G5TAPEA, [-45/45/0/90/- 30/30/0] s , 

Right Notch-Tip, at 61% S u lt- 
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da=0.20" 




90 degree ply 






da=0.20" 


-30 degree ply 




Figure A-2a - AS4\3501-6 Fiber Fracture and Delamination Illustration, 

Specimen G4TAPEA, [-45/45/0/90/-30/30/0] s . 

Both Notch-Tips, at 88% S u lt- 
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30 degree ply 


da=0.20" 



da=0.14" 


-30 degree ply 



da=0.29" 



90 degree ply 




0 degree ply 



da=0.25” 




Figure A-2b - AS4/3501 -6 Fiber Fracture and Delamination Illustration, 

Specimen G4TAPEA, [-45/45/0/90/-30/30/0] s . 

Both Notch-Tips, at 88% S u lt- 
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da=0.36" da=0.09" 



90 degree ply 




-30 degree ply 




da=0.16" 


30 degree ply 



da=0.09" 



da=0.32” 


0 degree ply 



da=0.18" 


Figure A-3a - AS4/938 Fiber Fracture and Delamination Illustration, 

Specimen D3AK5A, [-45/45/0/90/-30/30/0] s , 

Both Notch-Tips, at 90% S u it- 
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30 degree ply 


da=0.21 " 


da=0.16" 



da=0.26" 


-30 degree ply 



da=0.15" 



da=0.22" 


0 degree ply 



da=0.17" 



90 degree ply 




45 degree ply 




-45 degree ply 



45/-45 interface 


The data at this 
location was lost 
due to de-ply 
damage 


Figure A-3b - AS4/938 Fiber Fracture and Delamination Illustration, 

Specimen D3AK5A, [-45/45/0/90/- 30/30/0] s , 

Both Notch-Tips, at 90% S u it. 
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APPENDIX B 
FIBER BRIDGING 

Fiber Bridging Literature Survey 

Fiber bridging has been scrutinized and studied as a crack growth inhibitor since the 
1980's. Since then much research has been done to model and predict fiber bridging 
behavior. Initially, there were typically two types of fiber bridging models - the steady- 
state fiber bridging (SSFB) models and the generalized fiber bridging (GFB) models. 
Aveston et al. [B.l] and Budiansky et al. [B.2] used an energy balance approach in SSFB 
models to derive an expression for K m , the stress intensity factor (S.I.F.) for the matrix, in 
terms of composite microstructural parameters under conditions of steady-state cracking 
during monotonic loading. The steady-state S.I.F. is independent of crack length. A 
continuum fracture mechanics analysis is combined with a micromechanics analysis in a 
GFB model to derive S.I.F. solutions for matrix cracks of arbitrary size. The constraint 
due to the intact fibers in the wake of the matrix crack is idealized as an unknown closure 
pressure. The models developed by Marshall et. al., McCartney, and McMeeking and 
Evans [B.3-B.5] are all GFB models and are commonly referred to as the MCE, MC, and 
ME fiber bridging models, respectively. 

The works of Aveston, Budiansky, Marshall, etc., initiated and inspired further 
studies of fiber bridging. Sensmeier and Wright [B.6] studied the effects of fiber bridging 
on fatigue crack growth in titanium matrix composites. The framework for this analysis 
was the MCE [B.3] fiber bridging model. Further utilization of the GFB models by 
Bukuckas and Johnson [B.7] was in a study of matrix fatigue crack growth behavior in 
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center-notched titanium matrix composites. They assumed the intact fibers in the wake of 
the crack are idealized as a crack closure pressure in their use of the MCE, MC, and ME 
[B.3-B.5] fiber bridging models. Finally, Chan [B.8] presents a theoretical analysis that 
examines the effects of cyclic degradation of interface on fiber bridging of fatigue cracks 
in metal matrix or intermetallic matrix composites. He calculated frictional stresses on 
individual fiber/matrix interfaces using crack tip micromechanics and the fiber bridging 
models based on the works of Marshall et. al. [B.3] as well as Hutchinson Jensen [B.9]. 

Fiber bridging models that are somewhat independent of the MCE, MC, and ME 
models have been under development. Bao and Song [B.10] derived crack bridging 
traction laws that are based on a fiber pull-out analysis coupled with three proposed 
fiber/matrix interface assumptions. Yin [B.l 1] introduced a fiber bridging model based 
on crack closure tractions and applies a superposition to the stress intensity factor 
solution of a center-cracked tension specimen to provide a modified stress intensity factor 
which includes the effects of fiber bridging. 

In light of the idea that fiber bridging increases fracture toughness and can be 
construed as a crack growth resistance mechanism, much research has occurred to 
correlate the effects of fiber bridging with crack growth resistance curves (R-curves). 

Suo et. al. [B.12] developed crack closure tractions in the form of spring laws (linear and 
non-linear are compared) which are inferred from experimental delamination R-curves. 
Miyajima and Sakai [B.l 3] used the experimental R-curve to study fiber bridging where 
the fiber bridging tractions are estimated by the Dugdale approach. Similarly, Sakai etal. 
[B. 14] used the Dugdale approach to estimate fiber bridging tractions. Fiber pull-out and 
bridging processes in the wake of the propagating crack tip are discussed in relation to 
experimental R-curves. 

An experimental investigation of the role of fiber bridging in the delamination 
resistance was conducted by Spearing and Evans [B. 15]. The results were compared with 
fiber bridging models utilizing a softening traction law. This led to schemes for 
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predicting trends in delamination resistance with specimen geometry and crack length. 

Hu and Mai [B.16] uses a crack bridging theory which considers the difference in 
experimental and theoretical compliances to determine the fiber bridging stresses in the 
form of a crack closure or softening law. They showed that the delamination R-curve is 
consistent with the observation of fiber bridging in the delaminated region. Finally, 
influence of the bridging zone length on the resistance curve behavior was examined by 
Zok and Horn [B.17], Experiments are correlated with fiber bridging models and 
compared with R-curves. They demonstrated, with a model utilizing a crack closure 
pressure, that resistance curves for composites depend on both the absolute length of the 
bridging zone and the length of the bridging zone relative to the total crack length and 
specimen width. 

Fiber Bridging and Translaminate Fracture 

There are currently three types of fiber bridging; elastic fiber bridging, frictional fiber 
bridging, and pull-out fiber bridging. Figure B-l. Elastic fiber bridging is the case where 
the crack circumvents the fiber such that the fiber and matrix interface remain intact The 
interfacial shear strength, in frictional fiber bridging, is exceeded causing interfacial 
debonding and frictional stretching without fiber fracture. Finally, in pull-out fiber 
bridging, the fiber is shorter than the debond length and/or fractures within the debond 
length. For all three types of fiber bridging, previous works have stated and shown that 
fiber bridging is a crack growth resistance mechanism (toughening mechanism). 

However, all three types of fiber bridging are matrix cracks being bridged by intact Fibers. 

Consider the shear lag model [B.18, B. 19]. In the concept of shear lag, there is a 
region where interfacial shear stresses exceed the strength of the interface. It is in this 
region where there is relative sliding between the fiber and matrix. This often results in 
fiber bridging - fibers within the wake of the crack remain intact. Figure B-2. It can be 
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said then that fiber bridging is a toughening mechanism and is therefore a contributing 
factor in the R-curve behavior. The controlling mechanism for the case of shear lag is the 
high interfacial shear stress. 

The work herein focuses on translaminate fracture, not just matrix cracking alone. 

The bridging effects in translaminate fracture differ from the current fiber bridging 
configurations as is illustrated in Figure B-3. Typically, for translaminate fracture, there 
exists a fractured ply bridged by neighboring intact plies; usually a 0 degree ply bridged 
by off-axis plies. The ply bridging affects the redistribution of load into the neighboring 
plies. For a better understanding of this, it is a good idea to compare the behavior of 
metals to fiber reinforced composites. Figure B-4. We know that there are a lot of 
similarities in the mechanical behavior of metals to fiber-reinforced composites but their 
toughening mechanisms are different. For instance, the Dugdale plastic zone considers 
an effective crack longer than the physical crack. Crack edges in front of the physical 
crack carry the yield stress, tending to close the crack. The size of Aa is chosen such that 
the stress singularity is eliminated. It is proposed here that a similar phenomenon to the 
crack closure in the plastic zone is occurring in fiber reinforced composites, Figure B-5. 
Here, the fracture toughness is equal to the toughness found from the applied loading plus 
the additional toughness due to the bridging effects. The difference from the Dugdale 
approach is that the toughening mechanism and the calculations of the bridging effects 
are entirely different. Recent works by Poe [B.20, B.21] state that the general fracture 
toughness parameter, Q c in Figure 33, is solely a material parameter. However, if the 
structural effects of fiber bridging are present, the general fracture toughness parameter 
will be affected by its presence. 
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Figure B-l - Three (Current) Types of Fiber Bridging. 
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Current Fiber Bridging 
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Figure B-3 - Current Fiber Bridging Configurations and Translaminate Fiber Bridging. 
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Figure B-4 - Fracture Toughness Solutions for Metals. 
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Figure B-5 - General Fracture Toughness and Translaminate Bridging. 



APPENDIX C 

OPEN-HOLE TENSION MESH REFINEMENT STUDY 


Previous work [34] demonstrated the ability of the progressive damage model to 
accurately predict stiffness loss of open-hole IM7/5260 composite laminates loaded in 
tension-tension fatigue. Figure C-l. As part of the development of the residual strength 
methodology, a mesh refinement study of the open-holestrength cases was conducted. 

Four quarter-panel meshes were generated for the one inch wide and eight inch long 
open-hole specimens. Only about a third of the length is shown for illustrative purposes 
in Figure C-2. The result of this study is given in Figure C-3 which leads to the 
conclusion that the solution converges very well for the open-hole tension tests. Residual 
strength predictions for the open-hole laminates is given in Figure C-4 for completeness. 
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Mode I Matrix Cracking 
in the 90 Degree Plies 


Matrix Cracking in the 
+/- 45 Degree Plies 


Delamination 


Matrix Cracking 
in the 0 Degree Plies 
(Axial Splitting) 


Figure C-l - Tension-Tension Fatigue Damage in a Notched 
[0/45/-45/90] s IM7/5260 Laminate. 
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Mesh 1. 387 nodes 336 elements 



Mesh 2. 1071 nodes 992 elements 



Mesh 3. 2225 nodes 2112 elements 



Mesh 4. 2813 nodes 2688 elements 



Figure C-2 - Finite Element Meshes Used in the Convergence Study. 
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Figure C-3 - Mesh Refinement Study for the Residual Strength Predictions 
of the [0/45/-45/90] Laminate Open-Hole Geometry. 
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Figure C-4 - Predictions of Residual Strength. 
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APPENDIX D 

R-CURVES / CHARACTERISTIC CRACK LENGTHS 


The R-curves in Figures 34-45 were developed from the calculated effective crack 
growth, Aa. Another way to determine the effective crack growth is to simply measure it 
from the x-ray radiograph. This method is subjective because one needs to determine 
from the x-ray radiograph exactly what represents Aa. This may vary from person to 
person and from one x-ray radiograph to the next. An example is given to illustrate the 
procedure used to generate the R-curves, both for the measured Aa as well as for the 
closed form solution Aa. This example uses data from one of the experiments and is 
representative of all of the calculations involved in generating the R-curves in this study, 
including the predicted R-curves. 

Discontinuities along the load/cod plot in Figure D-l are labled A, B,C, and D. It is 
at these discontinuities that the specimen was unloaded and x-ray radiographs were taken 
just like the ones in Figures 14-16. The damage at the notch-tip was measured 
transversely from the tip of the notch toward the outer edge of the specimen with a 
machinists scale with divisions of 1/100 of an inch. The x-ray radiographs were enlarged 
to make this task easier. Table D-l illustrates the steps taken to arrive at the Aa 
measurements given in Figure D-l. 

To determine Aa from the closed form solution shown in Figure D-2, a straight line is 
drawn up the initial slope to aid in obtaining the initial point of nonlinearity. The load at 
this initial point is labeled Pj and is about 12.3 kips for this specimen. Given that 
E x =8.669 Msi, and the plot shows Sj=34,442 psi and codi=0.00492 in., then [cod*E x /4S]i 
has a value of 0.3096. Table D-2 demonstrates the final steps to calculating Aa and 
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fracture toughness. The failure values in Table D-2 were calculated for the point of 
catastrophic failure. These values are not represented by data points in the plots, they 
were only calculated for the fracture toughness values in Tables 5-8. The R-curves using 
the measured crack growth. Figures D-3 to D-l 1, include experimental and predicted 
plots. 


Table D-l Measured Characteristic Crack Growth Results 


Discontinuity 

X-ray 

Radiograph 

Magnification 

Magnified 
Aa (in.) 

Characteristic 

Crack 

Growth, 

Aa (in.) 

Stress, 
o (psi) 

Qc^tuf 

( VhT) 

A 

3.72 x 

0.40 

0.1075 

34,442 

0.1852 

B 

3.72x 

0.57 

0.1532 

40,808 

0.2343 

C 

3.72x 

0.97 

0.2608 

51,048 

0.3357 

D 

3.72x 

1.19 

0.3199 

54,229 

0.3810 

Ta 

ble D-2 Calculated Characteristic Crac 

k Growth Results 

Discontinuity 

COD (in.) 

Load, P (lb) 

Stress, 
G (psi) 

Characteristic 
Crack 
Growth, 
Aa (in.) 

Qc^tuf 

( VhT) 

B 

0.0062 

14,114 

40,808 

0.0163 

0.1875 

C 

0.0085 

17,656 

51,048 

0.0407 

0.2456 

D 

0.0099 

18,756 

54,229 

0.0679 

0.2736 

Failure 

0.0108 

20,261 

58,579 

0.0726 

0.2975 
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Figure D-l - Load/COD Plot for the Illustration of Crack Growth Measurements. 
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Figure D-2 - Load/COD Plot for the Illustration of Crack Growth Calculations. 
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Crack Growth, Aa (in.) 

Figure D-3 - Experimental and Model Generated R-Curves 
Using Aa Measured from X-Ray Radiographs. 
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Crack Growth, Aa (in.) 

Figure D-4 - Experimental and Model Generated R-Curves 
Using Aa Measured from X-Ray Radiographs. 
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Crack Growth, Aa (in.) 

Figure D-5 - Experimental and Model Generated R-Curves 
Using Aa Measured from X-Ray Radiographs. 
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Crack Growth, Aa (in.) 

Figure D-6 - Experimental and Model Generated R-Curves 
Using Aa Measured from X-Ray Radiographs. 
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Crack Growth, Aa (in.) 

Figure D-8 - Experimental and Model Generated R-Curves 
Using Aa Measured from X-Ray Radiographs. 
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131 


Crack Growth, Aa (in.) 

Figure D-9 - Experimental and Model Generated R-Curves 
Using Aa Measured from X-Ray Radiographs. 
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Crack Growth, Aa (in.) 

Figure D-10 - Experimental and Model Generated R-Curves 
Using Aa Measured from X-Ray Radiographs. 
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